Isospin dependence of relative yields of K + and K° mesons at 1.528A GeV 
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Results on K + and K° meson production in ||Ru + ||Ru and |q Zr + toZr collisions at a beam 
kinetic energy of 1.528A GeV, measured with the FOPI detector at GSI-Darmstadt, are investigated 
as a possible probe of isospin effects in high density nuclear matter. The measured double ratio 
(K + /K°)r u / (K + 1 K°)z r is compared to the predictions of a thermal model and a Relativistic Mean 
Field transport model using two different collision scenarios and under different assumptions on the 
stiffness of the symmetry energy. We find a good agreement with the thermal model prediction and 
the assumption of a soft symmetry energy for infinite nuclear matter while more realistic transport 
simulations of the collisions show a similar agreement with the data but also exhibit a reduced 
sensitivity to the symmetry term. 

PACS numbers: 25.75.-q, 25.75.Dw 



I. INTRODUCTION 



One of the main motivation to study relativistic heavy 
ion collisions at intermediate energies is to obtain infor- 
mation on the equation-of-state (EoS) for nuclear mat- 
ter at extreme conditions of pressure and density 
Within the high beam kinetic energy range of SIS (1-2 A 
GeV), the system could reach three times the normal nu- 
clear density and a temperature of about 60 MeV P, Q • 
Presently there are strong claims for a soft EoS based on 
the system size dependence of the K + yields 0, 0, [E @| 
and on the azimuthal dependence of the mean kinetic en- 



ergy of Z=l particles pj. However, when considering the 
excitation function of differential directed [1] and elliptic 
flow @ of charged particles in the incident energy range 
from 0.1 to 1.5 A GeV, the stiffness of the EoS appears 
less consistent, despite earlier preference for a soft EoS 
[l[ from similar measurements. 

The EoS of asymmetric nuclear matter expressed as 
the energy per nucleon as a function of baryonic den- 
sity pb and isospin asymmetry a is usually described as 
follows [HI [ID]: 

E(p B , a) = E(p B , a = 0) + E sym (p B )a 2 + 0(a 4 ), (la) 

l d 2 E( PB ,a) l N-Z 
E sym ( PB )=2 q^2 la=o, a ^J^Tz- (lb) 
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E S ym is the nuclear symmetry energy per nucleon and a 
the asymmetry parameter. 
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Predictions for the density dependence of E syr n based 
on various many-body theories diverge widely fiol [Til . 
EH, 0, Even the sign of the symmetry energy at 

baryonic densities above three times the normal nuclear 
density is still uncertain [l5| . The study of the isospin 
dependent part of the EoS is important to understand 
astrophysical processes like the structure of neutron stars 
and supernovae explosions jig, [lTf but also the structure 
of nuclei far off the valley of stability [H, [lj| . Recently, 
the momentum dependence of the symmetry potential 
was implemented in transport models in order to describe 
its effect in heavy-ion collisions more precisely [ly 1 , [2l| . 

Several observables, for example n/p and n~ /tt + yield 
ratios, are shown to be sensitive probes of the high den- 
sity behavior of the nuclear symmetry energy within 



131, [3]. Accord- 
H, the K+/K° 



transport model calculations 
ing to recent theoretical results 13, 
yield could be very sensitive to the isospin contribution 
to the EoS. For beam energies near kaon production 
threshold, kaon yields should carry information about 
the dense and hot system created during the collision 

ii i m, MM- 

In this work, we investigate the behavior of the sym- 
metry energy of the nuclear matter by comparing the 
relative production yields of K + and K° to theoretical 
predictions [13, [H, 127} • The study of the double ratio 
(K+ /K ) Ru /(K+/K )zr permits to isolate the isospm 
contribution on kaon production in systems (IfRu and 
40 Zr) which have the same atomic mass but a differ- 
ent number of neutrons. The FOPI collaboration has 
performed an experiment to maximize differences in the 
asymmetry (Eq.jTb]) with ^Ru + ^Ru and fgZr + ^Zr 
systems corresponding to asymmetry parameters a of 
0.083 and 0.167, respectively. The present results are ob- 
tained at 1.528A GeV beam kinetic energy (^/snn = 2.53 
GeV), slightly below the threshold energy of K + and K° 
production. Isospin effects measured at 400A MeV for 
the same colliding nuclei have been reported in [H, [2i| . 

The paper is structured as follows. Section II consists 
of a short description of the apparatus and the event char- 
acterization. The identification method of K + and the 
reconstruction method of K s are described in section III. 
In section IV the experimental K + /K° ratio measured 
in Ru and Zr systems is compared to theoretical model 
predictions. Finally, we summarize and give an outlook 
in section V. 



II. EXPERIMENTAL SET UP 

The experiment has been performed with the FOPI 
detector at the Heavy-Ion Synchrotron SIS of GSI- 
Darmstadt by using IfRu and IgZr beams of kinetic en- 
ergy of 1.528A GeV on enriched targets (>95%) of f 4 Ru 
and 40 Zr. The beam intensity was 3 x 10 4 ions/s, and 
the target thicknesses were 431 and 380 mg/cm 2 for Ru 
and Zr, which is equivalent to an interaction probability 
of 1%. More details on the experiment can be found in 



The FOPI detector is an azimuthally symmetric appa- 
ratus made of several sub-detectors which provide charge 
and mass determination over nearly the full solid an- 
gle. The central part of the detector, placed in a super- 
conducting solenoid, consists of a Central Drift Cham- 
ber (CDC) surrounded by plastic scintillators (Barrel). 
The forward part is composed of a wall of plastic scin- 
tillators (PLAWA and ZDD) and another drift chamber 
(Helitron) placed inside the super-conducting solenoid. 

For the present analysis, the PLAWA, the Barrel and 
the CDC were used. Particles measured in the CDC 
are identified by their mass, which is determined by the 
correlation between magnetic rigidity and specific energy 
loss. The Barrel and the PLAWA provide charge identifi- 
cation of the reaction products, combining time of flight 
(ToF) and specific energy loss information. More de- 
tails on the configuration and performances of the differ- 
ent components of the FOPI apparatus can be found in 

jaSmiHl. 

The events are selected by their centrality according 
to the charged particle multiplicities measured in the 
PLAWA and in the CDC. The results presented in the 
following are obtained for the most central 260 mb cross 
section, which corresponds to 7% of the total geometrical 
cross section. 



III. IDENTIFICATION AND 
RECONSTRUCTION OF KAONS 

A. Identification of K + 

K + mesons are identified by correlating their momen- 
tum determined with the CDC and their velocity mea- 
sured with the ToF Barrel. In addition, a cut on the 
laboratory momentum {pi a b < 0.5 GeV/c) is applied in 
order to reduce the contamination due to fast 7r + mesons 
and protons. More details on the K + identification are 
provided in [H, [34| . The phase space occupancy of K + 
in both systems (Ru + Ru and Zr + Zr) corresponds to 
a covering of the laboratory polar angle Ql from 39° to 
130°. The efficiency of K + identification is the same for 
the two systems and the fraction of measured kaons is 
about 12% of the total number produced in Air. 

An example of the K + mass spectrum is shown in 
Fig. [U for Zr + Zr collisions. The distribution was fit- 
ted with exponential and Gaussian functions in the mass 
region between 0.3 and 0.7 GeV/c 2 . The background is 
estimated with the exponential function while the signal 
is described by the Gaussian function. The signal af- 
ter background subtraction is plotted in the lower panel 
of Fig. [0 Within an interval of ±2cr, about 1180 K+ 
are identified for the applied set of cuts in the analy- 
sis. The signal-to-background ratio is about 2.2 and the 
mass width (a) of 32 MeV/c 2 is dominated by the de- 
tector resolution. The statistical error of the signal is 
evaluated with the covariant matrix by using the error 
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FIG. 1: Mass spectra of K + , the solid and dashed lines de- 
note the signal and background, respectively (upper panel). 
The lower panel shows the signal after background subtrac- 
tion. The following characteristics of the signal are reported: 
number of counts in the signal (S), significance (SIGNIF) and 
signal-to-background ratio (S/B) at ±2cr. The fitted param- 
eters of the Gaussian function as the mean value of the mass 
peak (MEAN) and the Gaussian width (a) are also mentioned. 



on Gaussian parameters. 

The systematic error of this measurement is deter- 
mined by applying different cuts on sensitive variables 
used to identify K + such as the matching quality be- 
tween the CDC and the Barrel. Also, another description 
of the background by using a polynomial function was 
used. After the determination of the number of counts 
in the mass spectra, a normalization with respect to the 
number of events was applied in order to directly com- 
pare the two colliding systems. Finally, the ratio of the 
production yields of K + measured in the two systems is: 

K t,J K Zr = L06 ± 0.07(stat.) ± 0.09(syst.) (2) 



B. Reconstruction of K° 

Due to K° — K° mixing, K^ and Kg are produced 
in equal amounts. As a consequence of their lifetime 
and decay channels, the K° can not be measured with 
the FOPI detector whereas the K s can be reconstructed 
via the weak decay into n~ and 7r + (branching ratio = 
68.6%, ct = 2.68 cm). The lifetime of Kg is long enough 
to permit the reconstruction of its vertex decay in the 
CDC. Therefore it is possible to significantly improve the 
signal-to-background ratio by rejecting particles coming 
from the primary vertex. The phase space occupancy of 
Kg in both systems (Ru + Ru and Zr + Zr) corresponds 
to a covering of the laboratory polar angle 0l from 32° 
to 140°. 



The Kg are reconstructed from identified pions by ap- 
plied a cut on the laboratory momentum (pi a b > 0.08 
GeV/c) in order to reject pions which are spiraling in 
the CDC. The distance of closest approach (DC A) be- 
tween tracks and the primary vertex in the transverse 
plane (DCA > 0.7 cm) is used to enhance the fraction 
of particles coming from a secondary vertex in the data 
sample. The efficiency of K s reconstruction is the same 
for the two systems. With the above mentioned cuts, we 
reconstruct 8% of all Kg produced. 

The invariant mass of the particle pair is calculated 
from the four-momenta of the pions at the intersection 
point. The reconstruction of K s measured for the Zr + 
Zr system is shown in Fig. 
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FIG. 2: Invariant mass spectra of tt^tt + pairs. The solid line 
and crosses denote the combinatorics and the scaled mixed- 
event background, respectively (upper panel). The lower 
panel shows the signal after background subtraction. 

The combinatorial background is obtained by using the 
event- mixing method (3o| . For that purpose, the two de- 
cay particles are taken from two different events which 
have the same centrality selection. In addition, the two 
events are rotated into the reaction plane in order to 
have the same reference system for both particles. This 
latter is estimated event by event according to the stan- 
dard transverse momentum procedure [36| . The mixed- 
event background is normalized in order to fit the com- 
binatorics in the invariant mass range 0.35-0.45 GeV/c 2 . 
The shape of the resulting mixed-event background de- 
scribes the combinatorial background and is indicated by 
crosses in Fig. [2] (upper panel). The vertical bars of the 
crosses correspond to the statistical errors. After back- 
ground subtraction (Fig. ® lower panel), the remaining 
peak in the mass spectra is fitted with a Gaussian func- 
tion. Within an interval of ±2cr, about 320 K s are found 
for the applied set of cuts in the analysis. The signal-to- 
background ratio is about 1.1 and the width (a = 11.4 
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MeV/c 2 ) of the distribution reflects the detector resolu- 
tion. 

The statistical error is determined by using the same 
method as for K + . In order to estimate the systematic 
error of this measurement we apply different sets of cuts 
by varying the selection criteria of the relevant quantities 
discussed in the beginning of this section. Finally, we 
obtain the following ratio of the production yields of K° 
for Ru + Ru and Zr + Zr systems: 

k rJ k Zt = °- 94 ± 0.12(stat.) ± 0.06(syst.). (3) 



IV. MODEL COMPARISON 

The experimental double ratio obtained in this analy- 
sis is: 



K Ru/ K Zr (K + /K°) Zr - 



(4a) 



(K+/K°) Ru 
(K+/K°) Zr 



1.13 ± 0.16(stat.) ± 0.12(syst.). (4b) 



Note that by using this double ratio, the effect of differ- 
ent phase space occupancies for K° and K + mesons are 
cancelled out (Eq. Haj) . 




DATA THERM. NL NLp NLp5 

Data Vs Models 

FIG. 3: Experimental ratio (K + /K°) Ru /(K + /K°) Zr (star) 
and theoretical predictions of the thermal model (cross) and 
the transport model with 3 different assumptions on the sym- 
metry energy: NL (circles), NLp (squares) and NLp<5 (trian- 
gles) . The INM and HIC calculations are represented by open 
and full symbols, respectively (see text for more details). The 
statistic and systematic errors are represented by vertical bars 
and brackets, respectively. 

The experimental result is compared to a thermal 
model prediction [2(| (left part of Fig. [3]) using a temper- 
ature of T = 52 MeV and a baryonic chemical potential 



of (iB = 756 MeV. It is important to stress that the ther- 
mal model assumes full stopping and mixing of the two 
nuclei and ensures isospin conservation via a chemical po- 
tential (jj, Ja = -7.1 MeV and -13.5 MeV for Ru + Ru and 
Zr + Zr systems, respectively). Also included in Fig. [3] 
(right part) are two sets of calculations based on Rela- 
tivistic Mean Field theory (RMF) [H [3, H3, H3 . One 
assumes infinite nuclear matter (INM, open symbols) us- 
ing a temperature of T = 60 MeV and a baryonic density 
of pb — 2.5po ,22] and the second describes a finite sys- 
tem as realized in heavy ion collisions (HIC, full symbols) 
[27l ]. Three different assumptions for the symmetry en- 
ergy are used in the transport model calculations. The 
reference ratio obtained from a non-linear Walecka model 
[l4| (NL, illustrated by circles in Fig. [3]) represents the 
prediction with a soft symmetry energy, only the kinetic 
contribution is taken into account. Assumptions labelled 
NLp (squares) and NLpS (triangles) show the expecta- 
tion of the model including the vector part (p) and both 
vector part plus scalar part (p + S) of the isovector mean 
field, which correspond to a gradually stiffer symmetry 
energy [ID, 0, H2| • The corresponding value of the sym- 
metry energy is reported on Fig. [3] (right part) for the 
INM case. HIC calculations have been performed explic- 

Zr systems, while INM 



itly for URu + ^Ru and |gZr - 
calculations are taken from \21 



The prediction of the thermal model is in perfect agree- 
ment with the data. It is worth to point out that the ther- 
mal model calculations reproduce very well the measured 
values of the single ratios of K + and K° in Ru and Zr sys- 
tems (Eq. [5] and Eq. [3J respectively). In the INM calcu- 
lations the system is in equilibrium (chemical, thermody- 
namical, isospin) and with a soft symmetry energy (NL, 
no isospin dependent collision term) the model prediction 
can be directly compared to the thermal model calcula- 
tion. As expected, both predictions are in agreement. 
Furthermore, the transport model permits to vary the 
stiffness of the symmetry energy as it is shown in Fig. [3J 
(right part). An increase of the (K+ /K°) Ru /(K+ /K°) Zr 
ratio as a function of the symmetry energy is observed for 
both the INM and the HIC calculations. With the INM 
scenario, a large enhancement (~70%) is predicted when 
going from a soft (NL) to a stiff (NL pS) EoS. The com- 
parison with the experimental point shows that our result 
is in favor of a soft symmetry energy. The isospin de- 
pendence, however, is dramatically reduced, when treat- 
ing the EoS within transport model calculations. Taking 
into account the experimental errors, the measurement 
does not permit to distinguish between the three differ- 
ent assumptions introduced in the HIC scenario. The 
sensitivity on isospin dependence of this observable is re- 
duced to 5% between the two most extreme cases (NL 
and NLpS). In contrast to the static INM calculations, 
in the HIC scenario two important dynamical effects oc- 
cur: fast neutron emission (mean field effect) and trans- 
formation of neutron into proton in inelastic channels 
(no-chemical equilibrium). This appears to be the ma- 
jor reason for the reduction of the isospin effect on kaons 
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ratios. Finally, for kaons measured at threshold energy 
from systems with asymmetry parameters differing by a 
factor two, the isovector contribution in nucleon-nucleon 
interaction influences very marginally the relative yields 
of K + and K°. Recent theoretical calculations 37[ sug- 
gest that a measurement of the double ratio at lower 
energies and with systems exhibiting a large difference 
in their asymmetry might increase the sensitivity on the 
isospin dependence of kaon production. On the experi- 
mental side, measuring a sample of kaons with sufficiently 
high statistics at sub-threshold energies will be a consid- 
erable challenge. 



CONCLUSIONS 



with the predictions of the thermal model and of the 
RMF model for infinite nuclear matter with a soft sym- 
metry energy. The RMF model for a collision system 
(HIC) shows a similar agreement with the data but the 
sensitivity to the symmetry term decreases by more than 
one order of magnitude. Within the model described 
above, the isovector contribution of nucleon-nucleon in- 
teraction influences weakly the relative yields of kaons 
measured at threshold energy and for systems differing 
in the asymmetry term by a factor two. A good under- 
standing of this phenomenon together with the present 
measurement will have a strong impact when designing 
new experiments aiming at the isospin dependence of the 
nuclear equation-of-state. 



We have reported the first measurement of the relative 
production yield of K + and K° mesons in Ru + Ru and 
Zr + Zr systems at a beam kinetic energy of 1.528A GeV 
which is close to the kaon production threshold. 

The experimental (K+/K ) Ru /(K + /K ) Zr ratio is 
compared to a thermal model and a model based on the 
Relativistic Mean Field theory using three different as- 
sumptions on the symmetry energy and two different col- 
lision scenarios. The measurement is in good agreement 
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